Electronic structure and the glass transition in pnictide and chalcogenide 
semiconductor alloys. Part II: The intrinsic electronic midgap states 
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We propose a structural model that treats in a unified fashion both the atomic motions and 
electronic excitations in quenched melts of pnictide and chalcogenide semiconductors. In Part I 
(submitted to J. Chem. Phys.), we argued these quenched melts represent aperiodic ppcr-networks 
that are highly stable and, at the same time, structurally degenerate. These networks are character- 
ized by a continuous range of coordination. Here we present a systematic way to classify these types 
of coordination in terms of discrete coordination defects in a parent structure defined on a simple 
cubic lattice. We identify the lowest energy coordination defects with the intrinsic midgap electronic 
states in semiconductor glasses, which were argued earlier to cause many of the unique optoelectronic 
anomalies in these materials. In addition, these coordination defects are mobile and correspond to 
the transition state configurations during the activated transport above the glass transition. The 
presence of the coordination defects may account for the puzzling discrepancy between the kinetic 
and thermodynamic fragility in chalcogenides. Finally, the proposed model recovers as limiting 
cases several popular types of bonding patterns proposed earlier, including: valence-alternation 
pairs, hypervalent configurations, and homopolar bonds in heteropolar compounds. 



I. INTRODUCTION 

In the preceding article,— we presented a chemical 
bonding theory for an important class of pnictogen- and 
chalcogen-containing quenched melts and glasses. These 
materials exhibit many unique electronic and optical 
anomalies^ not found in crystals, and also are of great 
interest in applications such as information storage and 
processings^ We argued these materials can be thought 
of as aperiodic ppc-networks made of deformed-linear 
chains that intersect at atomic sites at nearly right an- 
gles. Extended portions of the chains exhibit a perfect al- 
ternation of covalent and weaker, secondary bonds, even 
though the lattice as a whole is aperiodic. This bond al- 
ternation results from a symmetry breaking of a parent, 
simple cubic structure with octahedral coordination; this 
structure is uniformly covalently bonded. It is the inti- 
mate relation of the amorphous lattice to its covalently 
bonded parent structure that allowed us to rationalize, 
for the first time to our knowledge, two seemingly contra- 
dicting features of a bulk glass, i.e., its relative stability 
and structural degeneracy^ 

Yet the argued presence of the degeneracy of the ppa- 
network does not, by itself, guarantee that the network 
can be realized as an equilibrated supercooled liquid or 
a quenched glass: For instance, elemental arsenic can be 
made into an amorphous film but does not vitrify readily. 
To be a liquid, the network should contain a large equilib- 
rium concentration of structural motifs corresponding to 
the transition states for activated transport in quenched 
melts. 

The goal of the present article is to identify the bonding 
patterns of the transition-state structural motifs at the 
molecular level and describe the rather peculiar midgap 
electronic states that are intrinsically associated with 
such motifs. It will turn out that these electronic states 



are responsible for many of the aforementioned electronic 
and optical anomalies of amorphous chalcogenide alloys. 

Identification of structural motifs in vitreous materials 
based on local coordination is difficult because the usual 
concept of coordination, which is not fully unambiguous 
even in periodic lattices, becomes even less compelling 
in aperiodic systems. Alternatively, one can try to clas- 
sify such local motifs in terms of deviation from a puta- 
tive reference structure, while assigning a corresponding 
energy cost. Such deviations could be called "defected 
configurations." However, in the absence of long-range 
order, denning a reference structure in vitreous systems 
is, again, ambiguous. To make an informal analogy, is 
there a way to identify a typo in a table of random num- 
bers? The answer, of course, depends on the presence 
and specific type of correlation between the random num- 
bers. The random first order transition (RFOT) theory^ 
dictates that glasses do form subject to strict statistical 
rules prescribed by the precise degree of structural degen- 
eracy of the lattice, implying that correlation functions 
of sufficiently high order should reveal defects, if any. 
In fact, already four-point correlation functions in space 
capture the length scale of the dynamic heterogeneity in 
quenched melts Yet the only type of order in glasses 
that appears to be unambiguously accessible to linear 
spectroscopy is the very shortest-range order: The very 
first coordination layer is usually straightforward to iden- 
tify by diffraction experiments, while the strong covalent 
bonding between nearest neighbors is identifiable via the 
independent knowledge of the covalent radii of the perti- 
nent elements. Already the next-nearest neighbor bonds 
appear to exhibit a continuous range of strength and mu- 
tual angular orientation. 

We have argued* aperiodic pper-networks naturally ac- 
count for this flexibility in bonding in semiconductor 
glasses, while retaining overall stability. In doing so, we 
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proposed a structural model, by which aperiodic ppa- 
networks can be thought of as distorted versions of much 
simpler parent structures defined on the simple cubic lat- 
tice. There is no ambiguity whatsoever with defining co- 
ordination on a simple cubic lattice, thus allowing one to 
classify unambiguously the parent structures. 

We will observe that the most important and essen- 
tially the sole type of defect in ppcr-bonded glasses is 
singly over- or under-coordinated atoms. These defects 
turn out to host peculiar midgap electronic states with 
the reversed charge-spin relation, i.e., chemically they re- 
semble free radicals. We will argue that these defects in 
fact correspond to the electronic states residing on the 
high-strain regions intrinsic to the activated transport in 
semiconductor glasses proposed earlier by us.— On the 
one hand, these electronic states lie very deep in the for- 
bidden gap. On the other hand, they are surprisingly ex- 
tended, calling into question the adequacy of ultralocal 
defect models. This large spatial extent reveals itself by 
redistribution of the malcoordination over a large num- 
ber of bonds, in a solitonic fashion, and derealization 
of the wave function of the associated electronic state. 
The extended coordination defects are mobile, consistent 
with the conclusions of our earlier semi-phenomenological 
analysis' 3 that the peculiar electronic states are hosted by 
high-strain regions that emerge during activated trans- 
port in quenched semiconductor melts. 

The article is organized as follows. Section HTl reviews 
the conclusions of the RFOT theory on the concentration 
and spatial characteristics of the transition state configu- 
rations for activated transport in supercooled liquids and 
frozen glasses, and a general mechanism for the emer- 
gence of associated electronic states.— In Section IIII| a 
systematic classification of coordination defects in parent 
structures is carried out. In Section llVl we demonstrate 
that malcoordination defects in parent structures become 
delocalized in the actual, relaxed structure and make a 
connection with our earlier, semi-phenomenological anal- 
ysis of the solitonic states £ When making this connec- 
tion, we perform several independent consistency checks 
that the degeneracy of aperiodic ppc-networks is indeed 
compatible with the degeneracy of actual semiconductor 
alloys. Lastly we argue that the electronic states make a 
temperature-independent contribution to the activation 
barrier for liquid reconfigurations, which helps explain 
the apparent disagreement between the thermodynamic 
and kinetic fragilities in chalcogenides. 

We have alluded to two vulnerabilities of ad hoc defect 
theories, i.e., the difficulty in defining a putative refer- 
ence structure and the presumption of defects' being ul- 
tralocal at the onset. Conversely, such ad hoc theories 
do not explain self-consistently how the defects combine 
to form an actual, quite stable lattice. The present ap- 
proach resolves this potential ambiguity. First, since the 
defects are defined on a specific lattice in the first place, 
the question of their coexistence in a 3D structure is au- 
tomatically answered. Second, as discussed in Section fVl 
we will see how several popular defect theories proposed 



much earlier on phenomenological grounds are naturally 
recovered as the ultralocal limit of the present picture. 

II. BRIEF REVIEW OF THE RFOT THEORY 
AND PREVIOUS WORK 

Below we outline the minimum set of notions from the 
Random First Order Transition (RFOT) theory of the 
glass transition and other previous work, as necessary 
for the subsequent developments. Detailed reviews of 
the RFOT theory can be found elsewhere^*^ 

Mass transport in liquids slows down with lowering 
temperature or increasing density because of increasingly 
more frequent molecular collisions. At viscosities of 10 
Poise and above, however, the transport is no longer dom- 
inated by collisions but, instead, is in a distinct dynami- 
cal regime: The equilibrium liquid density profile^ is no 
longer uniform, but instead consists of sharp disparate 
peaks^ whereby each atom vibrates around a fixed lo- 
cation in space for an extended period of time. In other 
words, the liquid is essentially an assembly of long-living 
structures that persist for times exceeding the typical re- 
laxation times of the vibrations by at least three orders 
of magnitude.— Under these circumstances, mass trans- 
port becomes activated: Atoms move cooperatively via 
barrier-crossing events whereby the current long-lived, 
low free energy aperiodic configuration transitions lo- 
cally to another long-lived, low free energy configuration. 
The multiplicity of alternative aperiodic configurations 
is quantified using the so called configurational entropy: 
A region of size N particles has e s " N ^ kB distinct struc- 
tural states, where s c is the configurational entropy per 
per rigid group of atoms, often called the "bead."— The 
RFOT theory predicts that the configurational entropy 
at the glass transition is ~ 0.8fcs per beadfil 

Individual atomic displacements during the transi- 
tions are small, i.e., about the vibrational displace- 
ment at the mechanical stability edge, which is often 
called the Lindemann length <1l^ Although the acti- 
vated transport regime is usually associated with super- 
cooled liquids ) 16 i 17 many covalently bonded substances, 
such as Si02, exhibit activated transport already above 
the melting point 

During an activated reconfiguration, two distinct low- 
energy structures must coexist locally, implying that a 
higher free-energy interface region must be present. Be- 
cause the interface separates aperiodic arrangements, it 
has no obvious structural signature; non-linear spec- 
troscopy is generally required to detect the interface. The 
transitions proceed in a nucleation-like fashion; they are 
driven by the multiplicity of the configurations, subject 
to the mismatch penalty at the interface yyN i 18 i 19 

F(N) = 7VA — Ts c N, (1) 

see Fig. [I] Here, TV gives the number of beads con- 
tained within the nucleus and the coefficient 7 — 
%V3Trk B T\n(a 2 /d 2 L ire)il2° The lengths a and d L ~ 0.1a 
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FIG. 1. Typical nucleation profile for structural reconfigu- 
ration in a supercooled liquid from Eq. flj and its surface 
energy component ^\/~N : , normalized by the typical barrier 
height, which reaches (35 — 37)fcsT at T — T g . 



region of size N* , consistent with the quench being a 
higher free energy state than the corresponding crys- 
tal. In fact, th e to tal mismatch penalty in a sample 
of size N is j\/l^*(N/N*) = Ts c N and thus equals 
the enthalpy that would be released if the fluid crystal- 
lized at this temperature, save a small ambiguity stem- 
ming from possible differences in the vibrational en- 
tropies. The cooperativity size N* grows with the de- 
creasing configurational entropy^ (see Eqs. ([3]) and ((U): 
N*(T) ~ N*{T g )[(T g - T K )/(T - T K )] 2 . Still, it reaches 
only a modest value of 200 or so at the glass transition on 
one hour time scalej 14 ' 20 N*{T g ) = 200 corresponds to a 
physical size £ = a(A r *) 1 / 3 ~ 6.0 a, i.e. almost universally 
about two-three nanometers. This important prediction 
of the RFOT theory has been confirmed by a number 
of distinct experimental techniques 2 ^— and recently, by 
direct imaging of the cooperative rearrangements on the 
surface of a metallic glasses Now, the resulting concen- 
tration of the domain walls near the glass transition is, 
approximately, 



denote the volumetric size of a bead and the Lindemann 
displacement^! respectively. The reconfiguration time 
corresponding to the the activation profile in Eq. ([1]): 
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20, 



(5) 



r = T e Ft ' kBT = e^' ATs ' 



(2) 



grows in a super- Arrhenius fashion with lowering the 
temperature because of the rapid decrease of the con- 
figurational entrop y 22 ' 23 



s c ~ Ac p T g (l/T K - 1/T), 



(3) 



where Ac p — Tds c / 8T\t is the heat capacity jump at 
T g , per bead. The quantity Tk, often called the Kauz- 
mann temperature or the "ideal glass" transition tem- 
perature, corresponds to the temperature at which the 
configurational entropy of an equilibrated liquid would 
presumably vanish: Hereby the log number of alterna- 
tive configurations would scale sub-linearly with the sys- 
tem size. While the temperature Tk appears in certain 
meanfield models* 2 ^ the ideal-glass state, if any^ would 
be impossible to reach in actual liquids because of the 
diverging relaxation times, see Eq. @. 
The nucleus size 



N* = ( 7 /Ts c ) 5 



(4) 



where F(N*) = 0, is special because it corresponds 
to a region size at which the system is guaranteed to 
find at least one typical liquid state. As a result, the 
size N* corresponds to the smallest region that can re- 
configure in equilibrium. Each transition requires work 
jy/N to create and grow the interface, at the expense 
of relaxing the old interfaces. Thereby, the total num- 
ber of interfaces remains constant on average. One 
thus concludes that the liquid harbors one interfacial re- 
gion with an associated excess free energy 7V N* per 



In the rest of the article we will assume the material 
is just above its glass transition temperature T g , so that 
it represents a very slow, but equilibrated liquid. Be- 
low T gi the lattice remains essentially what it was at T g , 
save some subtle changes stemming from the decreased 
vibrational amplitudes and aging. 

The present authors have argued^ that the interfaces 
may host special midgap electronic states, subject to a 
number of conditions. These conditions are satisfied in 
many amorphous chalcogenides and pnictides.i The num- 
ber of the intrinsic midgap states is about 2 per interface, 
implying a concentration of about 2/£ 3 , where £ is the co- 
operativity length from Eq. ([5]) . The argument in Ref£ is 
independent from the present considerations and is based 
on the Random First Order Transition (RFOT) theory of 
the glass transition,- These intrinsic states are centered 
on under- or over-coordinated atoms and are surprisingly 
extended for such deep midgap states. In addition, the 
states exhibit the reverse charge-spin relation. The ex- 
istence of the intrinsic states allows one to rationalize a 
number of optoelectronic anomalies in chalcogenides in a 
unified fashion.— The interface-based electronic states in 
glasses are quasi-one dimensional and are relatively ex- 
tended along that dimension^ This quasi-one dimension- 
ality stems from the structural reconfigurations them- 
selves being quasi-one dimensional: Activated transitions 
between typical configurations occur by a nearly unique 
sequence of single-atom moves that are nearby in space. 
In the simplest possible Hamiltonian that couples elec- 
tronic motions to the atomic moves, by which a transition 
between two metastable states takes place, only relative 
positions of the atoms are directly coupled to the elec- 
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tronic density matrix^ 

n ° l = [(- 1 ) ne n4 l , S c n, S ~ W+l(4, S c n+l, S + H - C -) : 

n,s 

(6) 

where c\ s (c n s ) creates (removes) an electron on site 
n with spin s. The on-site energy, which reflects local 
electronegativity, is denoted with (— l) n e„. The electron 
transfer integral fn,n+i between sites n and n+1 depends 
on the inter-site distance d n ^ n +x. The sites are centered 
on beads, not atoms. The beads are numbered with index 
fi, 1 < n < N*, in the order they would move during the 
reconfiguration. 



III. CLASSIFICATION OF COORDINATION 
DEFECTS 

The structural model proposed by us in Ref. 1 presents 
a systematic way to classify amorphous structures and 
the associated electronic-structure peculiarities in vitre- 
ous ppa- net works. Without claiming complete general- 
ity, we will consider the following specific type of par- 
ent structures, which are largely based on Burdett and 
MacLarnan's model of (the crystalline) black phospho- 
rous and rhombohedral arsenic^ In this model, each 
vertex on a simple cubic structure is connected to exactly 
three nearest neighbors, where only right angles are per- 
mitted between the links. Two specific periodic ways to 
place the links according to this prescription correspond 
to the crystals of black phosphorous and rhombohedral 
arsenic. We have seen that this model has a significantly 
broader applicability, if one allows for distinct atoms and 
also vacancies in the original cubic lattice. For instance, 
let us take the parent structure for black phosphorus 
(Fig. 5(b) of Refi), place pnictogens and chalcogens in 
the rock salt fashion (as in Fig. 6(b) or 7 of Refi), and 
then omit every third pnictogen in a particular fashion 
(as in Fig. 7 of Refi). This procedure yields both the 
crystal structure and the stoichiometry of the archetypal 
pnictogen-chalcogenide compound Pn2Cli3 ("Pn"= pnic- 
togen, "Ch"=chalcogen). In the crystal, each pnictogen 
and chalcogen are three- and two-coordinated, respec- 
tively, thus conforming to the octet ruled The actual 
structure of a crystal or glass is thus viewed as a result 
of the following multi-stage procedure^ (1) start with the 
simple cubic structure with all vertices linked to all six 
nearest neighbors; (2) break the links and place vacan- 
cies to satisfy the stoichiometry and the octet rule. The 
resulting, lowered-symmetry structure is called the "par- 
ent" structure. To generate the actual structure from the 
parent structure, (3) shift the atoms slightly toward the 
linked vertices; (4) turn on the interactions by associat- 
ing an electronic transfer integral to each bond; and (5) 
geometrically optimize the structure, subject to the re- 
pulsion between the ionic cores and the variations in the 
electronegativity, if any. 

In the above prescription, it is generally impossible to 



satisfy the stoichiometry and the octet rule at the same 
time except in the case of periodic crystals, as we will 
see explicitly in a moment. As a result, the resulting 
parent structure will generally exhibit a variety of de- 
fects in the form of under- and over-coordinated atoms, 
such as three-coordinated chalcogens. Similarly, incor- 
porating elements of group 14, such as germanium, into 
the lattice of strictly three-coordinated vertices would 
introduce coordination defects. The just noted pres- 
ence of four-coordinated vertices presents a convenient 
opportunity to point out that the hereby proposed struc- 
tural model, though not unique, is special in the con- 
text of ppa-bonded solids in distorted octahedral geom- 
etry. Suppose that, instead, we started from a model 
where each vertex is /our-coordinated. The correspond- 
ing parent structures, if periodic, could be unstable to- 
ward a tetrahedrally bonded solid with the /3-tin or re- 
lated structured As a result, the present discussion is 
likely limited to compounds where elements from group 
14 are in the minority, except when these elements come 
in pairs with chalcogens, such as in GeSe or Ge 2 Sb 2 Te5. 
Such pairs are isoelectronic with a pair of elemental pnic- 
togens, and so both constituents of the pair are three- 
coordinated. This is not unlike, for instance, how GaAs 
forms the diamond (zincblcnde) structure, in which each 
atom is four coordinated. 

The three-coordinated parent structures on a simple 
cubic lattice, with right angles between the bonds, have 
another very special property^ In any lattice satisfy- 
ing this constraint, whether periodic or not, each bond- 
containing line is a strictly alternating bond/gap pat- 
tern. Informally speaking, one may think of the lattice 
as made of linear chains, each of which consists of white 
and black segments of equal length in strict alternation; 
the junctions between adjacent segments correspond to 
the lattice vertices, while the white and black segments 
correspond to no-link and link respectively. This ob- 
servation allows one to easily estimate the number of 
all possible distinct parent structures. Since there are 
only two distinct ways to draw a perfect bond/gap al- 
ternating pattern, the total number of the parent struc- 
tures in a sample of size L\ x L 2 X L% with no vacancies 
is 2( L i i 2+LiL3+i2i 3 )/a 2 ; w here a is the lattice spacing. 
Note that in contrast to Ref. 30 , we count as inequivalent 
two structures that can be obtained from each other by 
rotation, because they will both contribute to the phase 
space. We immediately observe that the total number of 
distinct parent structures is sub-thermodynamic in that 
it scales exponentially with the surface of the sample, not 
its volume. As a result, such defect-free structures can 
not contribute significantly to the library of bulk aperi- 
odic states. The only exception to this statement is the 
hypothetical Kauzmann, or ideal-glass state, in which 
the configurational entropy would strictly vanish^ 3 - One 
thus concludes that an equilibrium liquid must contain 
not only perfectly bond-alternating configurations, but a 
finite number of defected configurations per unit length 
in each line, where a "defect" consists of two or more 
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FIG. 2. Generation of singly malcoordinated atoms by a "cut- 
and-shift" procedure, see text. The filled and empty circles 
do not necessarily imply chemically distinct species, but are 
used to emphasize that it is the bond pattern that is shifted, 
not the atoms themselves. 



bonds (gaps) in a row. 

Furthermore, parent structures of certain common sto- 
ichiometrics, such as Pn2Ch3, must host a thermody- 
namic number of vacancies^ Adding vacancies to defect- 
less aperiodic parent structures should generally lead to 
the creation of coordination defects. Indeed, to maintain 
the stoichiometry, the vacancies should be spaced on av- 
erage by the same distance. On the other hand, one 
can construct a bond-breaking pattern with an arbitrar- 
ily large period, which can be made arbitrarily different 
from the average spacing between the vacancies. 

As a consequence of the deviation from the strict al- 
ternation pattern along the individual linear chains, a 
thermodynamic number of atoms in a representative par- 
ent structure must be either under- or over-coordinated. 
Note that after geometric optimization, defects in the 
parent structure will be generally mitigated or even re- 
moved. For example, at least one of the covalent bonds 
around a four-coordinated pnictogen will generally elon- 
gate in the deformed structure, in the pnictogen's at- 
tempt to attain the favorable valency 3. As a result, 
the overcoordination will be spread among a larger num- 
ber of atoms. An example of defect annihilation is when 
an extra bond on an overcoordinated atom annihilates 
with a missing bond on a nearby undercoordinate ver- 
tex. For these reasons, a defect in the deformed lattice is 
not well-defined. Here, we analyze several important de- 
fect configurations in the parent structure, whereby the 
local coordination is entirely unambiguous. Before pro- 
ceeding, we should note that the discussion is not limited 
to atoms of valence three and two. The above discussion 
applies even if a fraction of the vertices in the defectless 
parent structure shoud be four-coordinated, as would be 
appropriate for elements of group 14. Estimating the 
total number of distinct parent structures is no longer 
straightforward, however it is clear the number could be 
only lowered compared to the strict Burdett-McLarnan 
case. 
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FIG. 3. Examples of malcoordination defects of differ- 
ent orders on a pnictogen site and their relation to single- 
malcoordination defects along individual directions. A sim- 
ilar diagram can be drawn for chalcogens, which are doubly 
coordinated in a defectless parent structure, the bond angle 
equal to 90°. 



Let us begin the discussion of defected parent struc- 
tures with an ideal parent lattice with one singly over- 
coordinated vertex. This defect corresponds to two links 
in a row in one of the three chains crossing the vertex 
in question. Formally, one can obtain such a defect from 
a perfect Burdett-MacLarnan structure, for example, by 
choosing a vertex, removing the adjacent gap-segment 
from one of the chains and shifting the rest of the bond 
pattern on that side toward the vertex, see Fig. [2j Note 
that owing to pp7r-interactions, such shift generally in- 
vokes an energy cost or gain in addition to the cost of 
overcoordinating the chosen vertex, which we will dis- 
cuss later. One may obtain a singly undercoordinated 
vertex analogously, see Fig. [5] It is straightforward to see 
that more complicated, multiple-malcoordination config- 
urations can always be presented as superpositions of the 
simple, single-malcoordination configurations described 
above, see Fig.[3J For instance, a doubly-overcoordinatcd 
atom simply corresponds to overcoordination on two in- 
tersecting chains, at the vertex in question. Another 
common type of defect could be obtained by simply re- 
placing a gap in an ideal parent structure by a link. This 
defect is equivalent to two singly-overcoordinated atoms 
along the same line next two each other. Coordination 
defects of higher yet order can be constructed similarly. 
Note that bond angles other than 90° also amount to a 
coordination defect, even if the total number of bonds 
obeys the octet rule, see Fig. |3](e) . 

Single-malcoordination configurations can move along 
linear chains or make turns by bond switching, see Fig. 21 
whereby the defect moves by two bond lengths at a time. 
Direct inspection shows, however, that an attempt to 
make a turn at a pnictogen in Fig. [4] (such as atom 4) 
will create a double defect, such as in Fig. [31(e) , and is 
therefore likely to be energetically unfavorable. Gener- 
ally, only turns on atoms of the same parity as the defect 
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FIG. 4. Illustration of the motion of an overcoordination de- 
fect by bond-switching along a linear chain (from atom 6 in 
(a) to atom 4 in (b)) or making a turn (from atom 4 in (b) to 
atom 2 in (c)). 



center are allowed. Note that the distance between the 
nearest singly over- and under-coordinated atoms, along 
a chain, should be always an even number of bonds, if 
the malcoordination is along the very same chain. These 
"opposite" defects can be brought together and mutually 
annihilated. Conversely, one observes that two opposite 
defects separated by an odd number of bonds can not 
annihilate. A pair of such defects can be said to be topo- 
logically stable, see Fig. [5J In this Figure, we also illus- 
trate that this pair of defects is mobile, subject to the 
"turning" rule above, of course. In addition, the specific 
defect-pair configuration in Fig.[S]is special because it al- 
lows for defect motion without the creation of II-shaped 
bond motifs, in which three bonds constitute three sides 
of a rectangle. These motifs are somewhat energetically 
unfavorable, see the Appendix. Finally note that one 
of the nearest neighbors of the four-coordinated atom in 
Fig. [Sic) is not covalently bonded to any other atoms, im- 
plying the orientation of its bond with the central atom 
is relatively flexible. We will see below that owing to 
this extra flexibility and the said topological stability 
of the defect pair, the central atom is particularly un- 
stable toward sp 3 hybridization when positively charged. 
Nevertheless, since each of the individual defects can be 
removed, such sp 3 hybridized units do not represent a 
distinct type of defect. Conversely, the formation of sp 3 
hybridized atoms in ppcr-networks is reversible. 

We thus conclude that singly-malcoordinated atoms 
comprise the only distinct type of defect in the par- 
ent structure. The number of the motifs in the ac- 
tual deformed structure that originate from the coordi- 
nation defects in the parent structure is determined by 
the corresponding free energy cost. This cost consists 
of the energy cost proper of single-malcoordination de- 
fects, the energy of their interaction, and a favorable en- 
tropic contribution reflecting the multiplicity of distinct 
defect configurations. With regard to the interaction, 
it is easy to see that like defects should repel and de- 
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FIG. 5. A specific example of the motion of a topologically 
stable pair of an over- and under-coordinated defect, whereby 
the two defects move to the right while switching chains. 



fects of the opposite kind should attract. Indeed, the 
higher the over (under) coordination is, the larger the 
energy cost; conversely, mutual annihilation of a singly 
over- and under-coordinate defects is energetically favor- 
able. One may thus conclude that for both energetic 
and entropic reasons, there are no bound states of single- 
malcoordination defects to speak of that could be clas- 
sified as a distinct type of defect. The only such bound 
state is a topologically-stable pair of a singly over- and 
under-coordinated vertex, as in Fig. [5j 

We have pointed out earlier— that defectless parent 
structures can be formally represented as a subset of 
configurations of the 64-vertex model, which is the 3D 
generalization of the 6- vertex model of ice and 8-vertex 
model of anti-ferroelectrics. 3 — The presence of malcoordi- 
nation can still be implemented in the 64-vertex model, 
by assigning an appropriate energy cost to such coordi- 
nation. However, the presence of pp7r-interaction implies 
non-adjacent bonds interact directly. One concludes that 
parent structures of pertinence to ppc-networks are sub- 
ject to a more general class of models, such as the Ising 
model with next nearest interactions. 34 

In concluding the discussion of defect classification, we 
emphasize that stoichiometry-based vacancies in parent 
structures should be regarded not as defects, but as an 
intrinsic part of the parent structure. As in the exam- 
ple of the parent Pn2CIi3 structure^ such stoichiometry- 
based vacancies contribute to driving the distortion of 
the original simple cubic lattice and are minimized in the 
deformed structure, subject to competing interactions. 
Formally, this means the energy cost of local negative 
density fluctuations in the deformed structure that stem 
from the vacancies in the parent structure is a pertur- 
bation to the ppa bonding, similarly to the sp-mixing. 
Elementary estimates of the energy cost of a vacancy in 
a glass2£ show the cost is prohibitively high, dictating 
that an equilibrium melt above the glass transition will 
contain negligibly few voids of atomic size or larger. Sim- 
ilarly, homopolar bonds are not regarded as defects from 
the viewpoint of the present structural model. In fact, 
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FIG. 6. Central part of neutral (Asli2)2i chain, whose 
ground state contain a coordination defect and the associated 
midgap level. The calculations are performed on semiempir- 
ical level by MOPAC 36 with PM6 parametrization and the 
RHF/ROHF method, see Ref.— for full computational detail 
and Supplementary Material for comparison between the 
results of optimization for several chain lengths. 

we have argued^ a thermodynamic quantity of homopolar 
bonds in parent structures are required for the presence 
of transport in supercooled melts. These notions are con- 
sistent with the relatively low electronegativity variations 
in pper-bonded materials 1 and hence a weak applicability 
of Pauling's rules^ To avoid possible confusion, we re- 
iterate that the present argument strictly applies only to 
equilibrium melts or to glass obtained by quenching such 
melts. Amorphous films made by deposition generally do 
not correspond to an equilibrated structure at any tem- 
perature and, presumably, could host a greater variety of 
structural motifs. 



IV. THE INTRINSIC MIDGAP ELECTRONIC 
STATES 

To characterize the motifs arising in the deformed 
structure as a result of one singly undercoordinated atom 
in a parent structure, we consider first an isolated chain 
of ppa-bonded arsenic atoms passivated by hydrogens: 
(Asli2)„. For concreteness, the chain is oriented along 
the z axis. An infinite undistorted chain is Peierls- 
unstable. 1 The two ground state configurations of an in- 
finite (AsH2) n chain consist each of a perfect alterna- 
tion pattern of covalent and secondary bonds. A good 
approximation for these ground states can be obtained 
using a finite linear chain containing an even number of 
monomers, see the preceding articled 

Even if the chain is embedded in a lattice, each of these 
two ground states can still be regarded as the result of a 
Peierls transition, so long as the deviation from the strict 
octahedral coordination is weak. It is shown in Ref.- 
that the ppa-chahi portion of the full one-particle wave- 
function is an eigenfunction of the effective Hamiltonian: 

H = U vva + V f {E - W en v) _1 V (7) 

with the same eigenvalue E, where Ti. vva contains ex- 
clusively the p orbitals in question, T-L en v the rest of the 
orbitals, and V the transfer integrals between the two sets 
of orbitals. Using Eq. ([7]), the effect of the environment 
can be presented as an (energy-dependent) renormaliza- 
tion of the on-site energies and transfer integrals of the 
ppcr-network proper, see Appendix C of Refi The most 
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FIG. 7. Electronic energy levels of the defected chain 
(AsH2)2i chain from Fig. [6] full MO calculation (crosses) vs. 
one-orbital model with renormalized ppa integrals (circles). 
States below the gap are filled; the midgap state is half-filled. 

significant contributor to this perturbation is the com- 
peting sp interactions. The perturbation resulting from 
these competing interactions turns out to be sufficiently 
weak, as could be inferred from (a) its magnitude being 
close to the corresponding estimate using a perturbation 
series in the lowest non-trivial order; and (b) the bond 
morphology in the chain conforming to what is expected 
from a ppa- network. 1 

Now, in contrast to an even-numbered chain, the 
ground state of a (4rt + l)-membered open chain can be 
thought of as having one singly undercoordinated atom, 
see Fig.Hl while the ground state of a (4n + 3)-membered 
open chain contains one singly overcoordinated atom. Ei- 
ther of these "defected" chains is a neutral radical with 
an unpaired spin. The spectrum of the defected chain, 
shown with crosses in Fig. [7J exhibits a singly-populated 
state exactly in the middle of the forbidden gap. The 
electronic wave function of this midgap state, shown in 
Fig.0 is centered on the undercoordinated atom and ex- 
hibits a significant degree of derealization. Likewise, the 
deviation from the perfect bond alternation, although the 
strongest in the middle of the chain, is delocalized over 
several atoms. The partitioning of the electronic den- 
sity between the p z orbitals proper (73%), the arsenics' 
s-orbitals (21%), and the rest of the orbitals (6%) in- 
dicate that the interactions that compete with the ppa- 
bonding proper are significant but, nevertheless, may be 
regarded as a perturbation to the latter bonding, simi- 
larly to the perfectly dimerized chain>i Note that the po- 
sition of the nodes of the wave-function is easy to under- 
stand using elementary Hiickel considerations: Consider 
an odd-numbered chain of identical atoms, each hosting 
an odd number of orbitals. A moment thought shows 
that the resulting set of orbitals (in one electron approx- 
imation) has a half-filled non-bonding orbital, whereby 
the wave-function has nodes on every second atom. 

The derealization of the malcoordination is manifest 
in the broad sigmoidal profile of the ppa transfer integral 
as a function of the coordinate, see Fig. [§] Because the 
malcoordination is distributed among a large number of 
bonds in the deformed structure, the resulting structural 
signature of such a defect is far from obvious. In fact, 
the magnitude of the transfer integrals for the bonds ad- 
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site 

FIG. 8. The wave function squared of the midgap state of 
the (Asfi2)2i chain: the circles correspond to the arsenics' p z 
atomic orbitals (AO) (total contribution 73%), triangles As s- 
AO's (21%), crosses - the rest of AO's (6%). All contributions 
are nearly zero on every other atom; these very small values 
are omitted for clarity. The solid line corresponds to Eq. (|10fl . 

jacent to the overcoordinated atom in a (An + 3) chain 
is very similar to that for the undercoordinated atom 
in a (4n + 1) chain. This magnitude is approximately 
the average of the magnitude of the transfer integral for 
the covalent and secondary bond in a perfectly dimerized 
chain, see Fig. [S] The same comment applies to the cor- 
responding bond lengths in the first place, see Fig. II T t in 
view of the direct relationship between d and t. As a re- 
sult, the malcoordination is essentially fully "absorbed" 
by the chain. For instance, the distances between atoms 
1 and 19 in 19- and 21-member chains are 48.71 and 
48.70 A respectively. These notions emphasize yet an- 
other time the ambiguity of defining coordination or coor- 
dination defects in actual, deformed structures. Despite 
being relatively extended, the undercoordination defect 
could be thought of as separating two alternative states 
with perfect dimerization. This defect exhibits the re- 
verse charge-spin relation, because it is neutral when its 
spin is equal 1/2. Adding or removing an electron would 
yield a singly-charged defect with spin 0. These char- 
acteristics of the midgap state in the (AsH^)™ chain are 
entirely analogous to those of the midgap state in trans- 
polyacetylene, which is a classic example of a Peierls insu- 
lator. Until further notice, we will consider only neutral 
chains. Unlike the perfectly dimerized chain, a geometry 
optimized defected chain is somewhat curved on average; 
this effect is small, however. 

The surprisingly large spatial extent of the midgap 
state in the arsenic chain results from an interplay of 
several factors, similarly to how the spatial character- 
istics of the solitonic state in polyacetylene depends on 
the lattice stiffness, electronic interactions etc. Never- 
theless, it turns out that in both systems the extent of 
derealization is determined essentially by only two pa- 
rameters, namely the renormalized transfer integrals of 
the covalent and secondary bond in the perfectly dimer- 
ized chain, denoted with t and t' respectively. (In the 
case of polyacetylene, the transfer integrals are ppn, of 
course.) This simplification takes place because, similarly 
to the case of a perfectly dimerized chain, the effects of 



the competing intra-chain and external interactions on 
the midgap state largely reduce to the renormalization 
of the ppa interaction proper. To demonstrate this no- 
tion, we take the geometrically-optimized chain, extract 
only the ppa integrals, and renormalize them according 
to Eq. ([7|) while setting E equal exactly to the center of 
the forbidden gap. The spectrum of the resulting Hamil- 
tonian is shown with circles in Fig. [7] Setting E at the 
center of the gap should result in an error the greater, 
the further the state in question from the gap center. 
Yet one infers from Fig. [7] that the error in the spectrum 
is relatively small, consistent with the smallness of the 
perturbation. 

The spatial dependence of the renormalized ppa in- 
tegrals is shown in Fig. [§] with circles. In the con- 
tinuum limited 7 of the Su-Schrciffer-Heeger (SSH) model 
for midgap states in trans-polyacetylene,— the trans- 
fer integral depends on the coordinate according to the 
expression^: 

t + t' t — t' .„ , (n — nn\ ,„ N 
^ — + —(-1)^^], (8) 

where the soliton is centered on bond no and its half- 
width £ s is given by£ 

«.-|±|* m 

where d is the length of the bond in the parent structure. 
Using the renormalized values of the parameters t and t' 
from Refi, we obtain the solitonic profile shown by the 
solid line in Fig. [TJJ Note the difference between the two 
sets of transfer integrals is small. Finally, the electronic 
wavefunction in the same continuum limit is given by the 
function: 

4>l = (l/2)(d/6)cosh- 2 [(n - n )d/U (10) 

as shown with the solid lines in Fig. [5] Again, the agree- 
ment with the result of an electronic structure calculation 
using renormalized transfer integrals is good. The agree- 
ment is not expected to be perfect, however, because of 
the approximations inherent in the continuum limit^i of 
the SSH model. 

We thus observe that to describe both the ground state 
and defected configurations of linear chains from a ppa 
network, we can apply a Hamiltonian that contains ex- 
plicitly only the ppa interactions whose parameters are 
renormalized by other interactions. This resulting tight- 
binding Hamiltonian, which corresponds with the renor- 
malized Hamiltonian ([7|). reads as: 

Hel = [( -1 ) ne n4, s c n, S - *n,n+l(4 )S c «+l, s + H - C -) > 

(ii) 

where s (c n s ) creates (annihilates) an electron with 
spin s in the p z -orbital on atom n. The on-site en- 
ergy, which reflects the local electronegativity, is denoted 
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FIG. 9. Spatial dependence of the ppa transfer integral in the 
(AsEb^i chain. The actual transfer integrals from H ppa in 
Eq.Qare shown with crosses. The renormalized integrals from 
% in Eq. [7] are shown with circles. The solid line corresponds 
to Eq. ©. 



with (— l) n e„. The renormalized electron transfer inte- 
gral t njn +i between sites n and n+1 depends on the bond 
length dn t n+i. The latter is determined by the restoring 
force of the lattice, in addition to the Peierls symmetry 
breaking force stemming from the electronic degree of 
freedom in Eq. (fTTj) . The restoring force of the lattice 
includes both the intra- and extra-chain perturbations 
to the ppa interaction. We have seen in Refi that a 
portion of the intra-chain perturbation, namely the sp- 
mixing also contributes to the symmetry-breaking that 
results in bond-dimerization of the chain. Yet accord- 
ing to Eqs. ([!]) and (fTU]) . the participating competing 
interactions determine the spatial characteristics of the 
coordination defect largely through the values of only two 
parameters: t and t , in addition to the lattice spacing d. 

We are now ready to argue that the defect states cen- 
tered on singly malcoordinated atoms in pper-chains, as in 
Figs.JTElshould be identified with the intrinsic states pro- 
posed in Ref. 9 On the one hand, pper-bonded glasses meet 
the sufficient conditions for the presence of the intrinsic 
states, as argued in the preceding articled On the other 
hand, we have seen here a singly malcoordinated atom 
is the only type of defect present in ppc-bonded glasses, 
while all of its characteristics are identical to those estab- 
lished independently for the intrinsic states, including the 
reverse charge-spin relation and the relative derealiza- 
tion. In addition, both types of states emerge at the co- 
existence region of two distinct lowered-symmetry states, 
which originated from a higher symmetry state. In the 
rest of this Section, we discuss the microscopic insights 
arising from both the common aspects and complemen- 
tarity of the semi-phenomenological approach from Ref. 9 
(case 1) and the present, molecular model for the midgap 
states in ppcr-bonded glasses (case 2). 

Let us begin with the electronic Hamiltonian gov- 
erning the formation of the midgap states. In both 
cases, the electronic Hamiltonians are quasi-one dimen- 
sional, whereby the sites are situated along a deformed 
line in space. Although the Hamiltonians are identical 
notation-wise, they are distinct microscopically: In case 
1, Eq. ©, the sites refer to rigid molecular units, often 



called "beads." By construction, the beads are not signif- 
icantly disturbed by liquid motions and, conversely, inter- 
act only weakly with each other, comparably to the van 
der Waals coupling. Beads usually contain several atoms 
each; they are essentially chemically equivalent and fill 
out the space uniformly. The volumetric size of the bead 
is denoted with a*^ In contrast, in case 2, Eq. (ITT]) the 
sites are actual atomic orbitals which are generally not 
chemically equivalent. To check whether the two Hamil- 
tonians are consistent, we recall the RFOT theory's pre- 
diction that the smallest cooperativity size £ is just above 
2aji£ implying a supercooled melt can host at most one 
malcoordinated atom per ~ 4 beads at the onset of ac- 
tivated transport at a temperature T cr . The volumetric 
size a of the bead can be estimated based on the RFOT's 
prediction that the configurational entropy at the glass 
transition on one hour time scale is universally 0.8 per 
beadJ^ For As2Se3, this estimate yields about a half of 
the As2Se3 unit per bead, i.e. roughly one bead per ar- 
senic atom, and a ~ 4.0 Aj£ An inspection of the parent 
structure of As2Se3 in Fig. 8 of Ref. 1 demonstrates that 
it is possible to accommodate one malcoordination de- 
fect per 4 arsenics so that the defects are separated by at 
least two atoms. Since the defects are not immediately 
adjacent in space, the geometric optimization of the par- 
ent structure will not necessarily lead to their mutual 
annihilation. 

Next, we address the width of the interface and the 
related restriction on the spatial variation of the elec- 
tronic transfer integral. In Ref,-, we used the simplest 
spatial profile of bead displacement during a cooperative 
rearrangement to suggest that the soliton half-width £ s 
is bounded from below by £/2, which is about 3a at T g . 
Inspection of Fig. [5] indicates that the soliton, though 
still extended, is somewhat more narrow, i.e. £ s ~ 2a. 
(Recall also that d < a.) Using the condition 9 

t + t' 

^ < —a, (12) 

we then should revise the lower limit on the variation of 
the transfer matrix element to: 

t'/t > 0.3. (13) 

We have seen earlier— that the ppa bonded materials in 
question do indeed meet this revised condition very well. 
Finally we note that the prediction of Ref^ that the ex- 
tent of the wave function somewhat exceeds 2£ s is indeed 
borne out by the result in Fig. |SJ 

The RFOT theory makes specific predictions as to the 
free energy cost of a structural reconfiguration. For in- 
stance, the total mismatch penalty near the glass transi- 
tion, according to the discussion of Fig. [TJ is four times 
the typical barrier for the reconfiguration. At the glass 
transition, this corresponds to about 140fcsT g , or gener- 
ically about 3 to 4 eV. For the present microscopic pic- 
ture to be valid, this energy should exceed the energy 
cost of the structural defect, which includes the energy 
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FIG. 10. Central portion of a geometrically optimized double 
chain of hydrogen-passivated arsenics. The top and bottom 
chains, containing 19 and 17 arsenics, host an over- and under- 
coordinated arsenic respectively, c.f. Fig. [6] The numbers 
denote the bond lengths in Angstroms. 



of the associated lattice deformation and the electronic 
energy proper. We have established earlier— that the elec- 
tronic energy associated with the intrinsic state can be 
estimated via the optical gap E gap and the effective at- 
traction U e ff between like particles occupying the state: 
~ (-Egap — |£/ ff|) — 0.7i?g a p. This figure is about one quar- 
ter of the full interface energy predicted by the RFOT 
theory. Alternatively, according to the solution of the 
continuum limit of the SSH model, the cost of the intrin- 
sic state on an isolated defect can be expressed through 
its width £ s : 



yielding a similar figure to the above estimate. One ob- 
serves that, indeed, the present microscopic picture is 
internally consistent. There yet is another potential elec- 
tronic contribution to the defect energy. We have already 
mentioned that the creation or motion of a defect may by 
accompanied by a creation of II-like patterns of bonds in 
the parent structure; the latter are analyzed in the Ap- 
pendix. 

A chain hosting a singly under- and over-coordinated 
defect can always lower its energy via mutual annihi- 
lation of the defects, consistent with conclusions of the 
continuum approach of Rcf. 9 Two such defects on differ- 
ent chains, however, may be topologically stable against 
such annihilation, as in Fig. [5jc). Such a configuration 
presents a felicitous opportunity to examine the interac- 
tion of the defects in the deformed structure, for instance, 
by geometrically optimizing a double (4n + 1) — (An + 3) 
chain of hydrogen-passivated arsenics, see Fig. [TO] One 
observes that in the ground state of a resulting defect 
pair, the under- (over-) coordinated atom is negatively 
(positively) charged, so as to obey the octet rule, see also 
Section |Vj The resulting dipole moment of the chain is 
3.3 Debye; see the shapes of the HOMO/LUMO orbitals 
and electrostatic map of the chain in the Supplementary 
Material^ The As- As angles at the overcoordinated ar- 
senic are 110, 110, and 128 degrees. Note in perfect tetra- 
hedral coordination, these angles would be 109.5 degrees. 
In addition, the defects become more localized when in 
close proximity, as shown in Fig. 111! because of the par- 
tial cancellation of the malcoordination. Alternatively, 
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FIG. 11. Site dependence of the As-As bond length in sin- 
gle and double chains of hydrogen-passivated arsenics, as 
labeled by "1-chain" and "2-chain" respectively. The "1- 
chain, under" corresponds to a (AsH2)i7 chain similar to the 
(AsH2)2i chain from Fig. [6] and "1-chain, over" corresponds 
to a (AsH2)i9 chain. The double-chain entries pertain to the 
individual chains in Fig. 1101 "over" and "under" correspond 
to the top and bottom chain respectively. 



one may say the energy cost of this localization partially 
offsets the stabilization resulting from the binding of the 
defects. 

The extended nature of the coordination defects is con- 
sistent with their mobility: In the continuum limit, the 
pinning of a soliton is in fact strictly zero. Consistent 
with the latter notions, the energy gradients during ge- 
ometrical optimization of the (AsH2)„ chain above have 
very small numerical values. Notwithstanding this com- 
plication, an estimate for the barrier for soliton hopping 
event (by two bond lengths) can be made within the 
SSH model: Using the effective transfer integrals from 
Eqs. (J5J and ((9]) we estimate the energy of the transi- 
tion state configuration (corresponding to a move by one 
bond length) at approximately 0.1 meV, which is about 
10 _5 £ , gap . ppir interactions may also contribute to pin- 
ning, see the Appendix for quantitative estimates. While 
the above notions may apply relatively directly to neutral 
defects even in actual 3D lattices, the situation appears 
to be more complicated for charged solitons (which in fact 
greatly outnumber the neutral ones^) because of lattice 
polarization. This is work in progress. Finally, note the 
high mobility of the soliton seems particularly important 
in the context of cryogenic anomalies of glasses, which 
have been argued to arise from the low-barrier subset of 
the structural reconfigurations . 10 i 20 

As an additional dividend of the present discussion, 
let us see that the presence of an electronic - and hence 
T-independent - contribution to the RFOT's mismatch 
penalty between distinct low free energy aperiodic con- 
figurations allows one to partially explain the deviation 
of chalcogenides from the RFOT-predicted correlation 
between the thermodynamics and kinetics of the glass 
transition. Indeed, because of this T-independent con- 
tribution we expect corrections to the detailed tempera- 
ture dependence of the reconfiguration barrier T* from 
Eq. ([2]) and the value of the so called fragility coeffi- 
cient m = (l/T)91og(r)/9(l/r)|T=T s ^^ The RFOT 
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theory predicts that this coefficient should be equal to 
what one may call a "thcrmodynamically" determined 
fragility mthermo = 34.7 x Ac p (T g ) where Ac p (T g ) is the 
heat capacity jump at the glass transition per bead, 13 ^ 38 
Using the measured Ac p (T g ) per mole, and the bead 
size determination from the fusion enthalpy produces ex- 
cellent agreement for dozens of substances,— however 
among the outliers are the chalcogenides. For instance, 
for As2Se3, the measured m ~ 40, while the theoret- 
ically computed mthermo — 7.5. Determination of the 
bead size is ambiguous in As2Se3, because the corre- 
sponding crystal is highly anisotropic. In fact, in the 
fusion-entropy based estimate, the bead comes out to 
be a single atom, in conflict with the aforementioned 
requirement of chemical homogeneity and with Ref. 39 , 
where a bead was argued to contain at least two atoms. 
Using the RFOT's prediction that s c (T g ) ~ 0.8k B per 
bead, gives a better value mthermo — 17,— still not enough 
to account for the disagreement. Simple algebra shows 
that in the presence of a T-independent contribution to 
the surface tension, the estimate of the thermodynamic 
fragility from RefJ^ should be increased by a factor of 
[l + 2(SE)(T g -T K )/T K ], where ST is the relative contri- 
bution of the electronic states to the total surface energy 
at T g . Assuming the mentioned <5£ ~ 0.25 and using pa- 
rameters for As2Se3 from Ref.—, we get an increase by a 
factor of 1.45, which is a considerable correction in the 
right direction. Besides the experimental uncertainty in 
determination of Tk, we note also that more than one 
pair of solitonic states might accompany the transition. 
This would significantly increase the corresponding con- 
tribution to the surface tension and the fragility. 

We note that since (T g - T K )/T K cx 1/Ac p , this T- 
independent correction to the surface tension would be 
less significant for fragile substances. ((T g — Tk)/Tk J$ 1 
for known substances). This reflects the old notion 
that the temperature dependence of s c is the lead- 
ing contribution to the viscous slowdown in a super- 
cooled melt, except in strong liquids whose slowdown is 
nearly Arrhenius-like. Incidentally, the strongest liquids 
from the survey by Stevenson and Wolynes,— , namely 
GeC>2, and BeF2, and ZnCl2, do conform to the relation 
mthermo = 34.7 x Ac p (T g ). This is consistent with the 
present results, since we do not expect solitonic states in 
these wide-gap materials. 



V. THE CHARGE ON THE COORDINATION 
DEFECTS AND THE RELATION TO EARLIER 
DEFECT THEORIES 

We have argued previously- that the intrinsic midgap 
states should be typically charged: A half-filled defect is 
essentially a neutral molecule embedded in a dielectric 
medium and should be generally stabilized by adding 
a charge, because of ensuing polarization. To analyze 
distinct charge states of the coordination defects in the 
present approach, we will use the same methodology 
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FIG. 12. A compilation of the possible charged states 
of singly malcoordinated atoms. Here, "Ch" =chalcogen, 
"Pn"=pnictogen. Neutral states, not shown, imply dangling 
bonds and would be energetically costly. 



as in Sections |TTT] and |IV] above: We will consider dis- 
tinct charge states on the defects in the parent struc- 
ture, with the usual understanding that the correspond- 
ing motif in the actual, deformed structure will be signif- 
icantly delocalized. It is during the analysis of charged 
defects in the parent structure that we will be able to 
make a connection between the present approach and 
the much earlier specific proposals on the defect states 
in chalcogenides 4Mi 

The large spatial extent of the malcoordination-based 
states is a result of the lattice's attempt to mitigate the 
energy cost of what would have been a local defect in the 
parent structure. The converse is also true: According 
to Eq. qi4p. the energy cost of localizing malcoordina- 
tion would well exceed the energy density typical of an 
equilibrium melt. Thus in the analogies below between 
the present and earlier approaches, one should keep in 
mind that on the one hand, the earlier proposals ascrib- 
ing midgap states to specific local defects can be thought 
of as an ultralocal limit of the present theory. On the 
other hand, we see that this ultralocal limit is somewhat 
artificial in that it greatly overestimates the energy cost 
of the defects. 

Because there is only one basic type of the defect, 
namely a singly under- or over-coordinated atom, and 
three distinct charges: -1, 0, 1, the possible combina- 
tions of these variables are only few, see Fig. [T^l Already 
an elementary analysis yields that charged states are ex- 
pected to be stabilized, consistent with independent ar- 
guments from Ref. 9 Indeed, four of the resulting configu- 
rations, namely those listed in the "Lewis octet" sectors, 
satisfy the octet rules and thus are expected to be partic- 
ularly stable. Specifically, a positively charged overcoor- 
dinated chalcogen is chemically equivalent to a pnictogen 
and a negatively charged undercoordinated pnictogen is 
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equivalent to a chalcogen. These two configurations are 
expected to stabilize the distorted octahedral arrange- 
ment of the ambient lattice. On the other hand, the Pn|" 
configuration is unstable toward tetrahedral order, c.f. 
Fig. [5jc) and Fig. [lOl and as such, would tend to frus- 
trate the ppc-network. The situation with the entries in 
the "hypervalent" and "hypovalent" sectors appears less 
certain. First off, we have used the label "hypervalent" in 
the right bottom sector in reference to the formal electron 
count around the atom exceeding 8 and the low magni- 
tude of electronegativity variation in the compounds in 
question. The T-shaped bond geometry in the Ch^~ case 
was chosen by analogy with small molecules with the 
same electron count, such as CIF3. If the T-shaped bond 
geometry is indeed favored for this type of defect, there 
will be an additional penalty for a negatively charged 
defect to make a turn at a chalcogen. Now, the left up- 
per section represents a hypovalent configuration. In the 
worst case, these hypovalent configurations would favor 
a triangular bipyramidal arrangement, in which there is 
linear alignment for at least two bonds. We have not 
listed neutral defects in the table. These would amount 
to having dangling bonds and thus are deemed energeti- 
cally costly. Finally note that pairs of oppositely charged 
defects with complementary coordinations will either an- 
nihilate or be particularly stable, if their annihilation is 
topologically forbidden. This attraction is indicated by 
the double-ended arrows in Fig. [12J A specific realization 
of an attractive, topologically stable pair of defects in a 
deformed structure is shown in Fig. 1101 This configura- 
tion formally corresponds to a bound Pn^-Pn|" pair. 

Several specific configurations from Fig. [T^l includ- 
ing Cb-3, Chg , Ch+, and PnJ, have been proposed 
as relevant defect species earlier, based on the puta- 
tive presence of dangling bonds 4 -^i vector-type charge- 
density waves ^ and molecular dynamics studies using 
energies determined by tight-binding methods42. Specif- 
ically in the venerable approach of Kastner et al4^, 
pairs of defects, called valence-alternation pairs (VAP), 
could spontaneously arise in glass because of a pre- 
sumed instability of a pair of dangling bonds toward 
creation of two charged defects: 2Chg — ^Ch^+Ch^ or 
2Pn° — J-Pn^+Ch^. Despite common features with the 
VAP theory, the present picture is distinct in several 
key aspects, in addition to the aforementioned ones. We 
have seen in Section [TTT] the coordination defects are a 
necessary prerequisite for molecular transport and are 
themselves mobile. In other words, (the delocalized ver- 
sions of) the specific configurations from Fig. fTS] are sim- 
ply transient configurations arising during motions of the 
structural interfaces in a quenched melt or frozen glass. 



VI. CONCLUDING REMARKS 

We have proposed a structural model, by which the 
structure and electronic anomalies of an important class 
of vitreous alloys emerge in a self-consistent fashion. The 



model allows one to reconcile two seemingly conflicting 
characteristics of quenched melts and frozen glasses. On 
the one hand, these materials exhibit remarkable thermo- 
dynamic and mechanical stability, only slightly inferior 
to the corresponding crystals. On the other hand, these 
materials are also multiply degenerate thus allowing for 
molecular transport. The stability of the lattice stems 
from its close relationship with a highly symmetric, fully 
bonded structured The lattice's aperiodicity and the nec- 
essary presence of a thermodynamic number of transition 
states for structural reconfigurations dictate that the lat- 
tice exhibit a thermodynamic number of special struc- 
tural motifs, roughly one per several hundred atoms at 
T g . These motifs host midgap electronic states that are 
similar to solitonic states in trans-polyacetylene4^ The 
motifs correspond to coordination defects in the parent 
structure, which can be defined entirely unambigously. In 
contrast, in the actual materials, coordination is poorly 
defined because the interatomic distances and bond an- 
gles are continuously distributed. The presence and the 
characteristics of these defects have been established us- 
ing an exhaustive, systematic protocol, implying the clas- 
sification of the defects is complete. 

The actual molecular mechanism of both the stabil- 
ity of semiconductor glasses and the mobility of the de- 
fected configurations relies on the very special property 
of chalcogen- and pnictogen- containing semiconductors, 
namely their fully developed pper-bond network. Hereby, 
most atoms exhibit a distorted octahedral coordination. 
In each pair of collinear bonds on such an atom, one 
is fully covalent and the other is weaker, but still di- 
rectional. The two bonds are intrinsically related be- 
cause they exchange electron density, similarly to the 
donor-acceptor interactions. During defect transport, the 
stronger and weaker bond exchange identities, resulting 
in bond switching not unlike the Grotthuss mechanism of 
bond switching in water. The electronic states residing 
on these mobile coordination defects are thus an intrinsic 
feature of transport in ppc-bonded melts and glasses. In 
contrast, bond switching in known tetrahedrally bonded 
semiconductors does not occur because the bonds in these 
materials are homogeneously saturated throughout, ex- 
cluding the possibility of intrinsic states in these mate- 
rials, consistent with our earlier conclusions^ Note that 
all types of glassformers host the high strain interfacial 
regions that separate low free energy aperiodic config- 
urations and are intrinsic to activated transport. Yet 
glassformers exhibiting distorted octahedral coordination 
appear to be unique in their ability to host topological 
electronic states. The corresponding glasses are thus ex- 
pected to exhibit anomalies that we have associated with 
these electronic states^ including light-induced electron 
spin resonance. This expectation is consistent with the 
variation of the magnitude of light-induced anomalies 
across the Ge^Sei-a, series^ for 1/3 < x < 1/2, which 
display coordination ranging from tetrahedral (smaller x) 
to distorted octahedral (larger x). In this series, the octa- 
hedral ordering seems to correlate with the light induced 
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ESR and vice versa for the tetrahedral bonding^ 

The coordination defects exhibit topological features, 
such as stability against annihilation, if there is a mis- 
alignment in the motion of two defects. This peculiarity 
is related to another topological feature of these states, 
which is especially transparent in the continuum limit of 
Eq. ([H]): 9 H — —iva 3 d x + A(x)oi + e(x)ai. Here, at are 
the Pauli matrices, while —iv, e(x), and A(x) correspond, 
respectively, to the kinetic energy, local one-particle gap 
and variation in electronegativity. The local gap A(x) is 
space dependent and, in fact, switches sign at the defect, 
thus corresponding to a rotation of a vector (A, e). The 
orientation angle of this vector is the topological phase 
associated with the defect. The phases of defects trav- 
eling along different paths can not cancel, resulting in a 
stability against annihilation. 

Finally, what are the implications of the present re- 
sults for direct molecular modeling of amorphous chalco- 
genides and pnictides? Parent structures differing only 
by defect locations can be used as the initial configu- 
rations for MD simulations and may, conceivably, help 
find transition state configurations corresponding to re- 
alistic quenching rates. Yet even assuming the inter- 
actions can be estimated accurately and efficiently, the 
problem is still very computationally expensive, as high- 
lighted by the present work: The number of defectless 
parent structures, although sub-thermodynamic, is still 
very very large. We have seen it is plausible that these 
structures comprise the elusive Kauzmann state for these 
systems. 
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Appendix A: The pp7r-interaction 

When the coordination is exactly octahedral, n-shaped 
bonds patterns, such as those formed by atoms 2-1-6-7 
and 4-3-8-9 in Fig.fTBTb). give rise to an additional energy 
cost stemming from the pp7r-interaction, as explained in 
the following. It is partially owing to this cost, for in- 
stance, that the layers of AsSe3 pyramids in the actual 
structure of As2Se3 are so greatly puckered, see Figs. 7-9 
of Refi To elucidate the origin and estimate the strength 
of the pp7r-interaction, we consider for concreteness the 
configurations shown in Fig. 1131 periodically continued 
along the double chain. Only atomic p-orbitals lying 
in the plane of the figure are considered; these orbitals 
host two electrons per atom. The nonzero entries of the 
Hamiltonian matrix are indicated in Fig. I13f a) ; the nota- 



tions are from Refi The total energy of the configuration 
in Fig. [TST a) can be expanded in a series in terms of the 
ratios t' a /t a and t n /t a . In the second order in the expan- 
sion, the total energy is the sum of three contributions 
(assuming e = 0): The covalent bond energy is —2t a per 
bond, the secondary bond energy equals — (t' a ) 2 /2t a per 
bond, and the contribution due to the pp7r-interaction is 
given by 



per each vertical pair of atoms. One can similarly show 
that the energy of the configuration in Fig. [TBT b) has the 
same contributions as in (a) except the stabilizing con- 
tribution from Eq. (|A1I) . owing to the n shaped motifs. 
As a result, the configuration in Fig. IT37 b) is less en- 
ergetically favorable. According to Eq. (|A1[) . such a n 
motif generically costs about 0.1 eV; typical values of the 
transfer integrals can be found, e.g. in Table 1 of Rcf. 1 
This fact can be interpreted equally well either as stabi- 
lization due to the staircase motifs in configuration (a) or 
destabilization due to the II motifs in configuration (b). 

The presence of the pp7r-interaction generally affects 
the energetics of defect movement, as could be seen in 
Fig. [T31 The top chains in panels (a) and (b) can be 
obtained from each other by switching all the bonds, as 
could be achieved, for instance, by moving a malcoor- 
dination defect along the chain. As a consequence, the 
motion of such a defect is subject to a uniform potential. 
Now, each structural transition in a supercooled melt is 
accompanied by the creation of two defects of opposite 
charge and malcoordination*2 It follows from the above 
discussion that in the worst case scenario, the cost of 
separating two defects scales linearly with the distance, 
i.e. as TV 1 / 3 , where N is the size of the rearranged re- 
gion. The scaling of this cost is clearly subdominant to 
that of the mismatch penalty term ^N 1 / 2 from Eq. (JTJ) . 
Furthermore, at large enough N, it will be energetically 
preferable to remove the II patterns separating the origi- 
nal pair of defects by inserting another pair of malcoordi- 
nation defects between them. By this mechanism, which 
is not unlike the mechanism of quark confinement, the 
cost of the II patterns is limited by a fixed number which 
is comparable, but smaller than the typical energy of the 
interface, i.e., 7Y 1 / 2 . We suggest the actual effect of 
the Il-configurations is yet smaller. Indeed, distinct liq- 
uid configurations should have, on the average, the same 
number of Il-configurations, implying the defect move- 
ments are subject to a zero bias, on average. As a result, 
the magnitude of the overall bias for defect separation 
scales at most with a square root of the distance, leading 
to a correction to the energy cost of the reconfiguration 
that is proportional to (TV 1 / 3 ) 1 / 2 = iV 1 / 6 . This correc- 
tion is clearly inferior to the 7./V 1 / 2 term. 
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(a) 



(b) 



o- 



FIG. 13. Double chain interacting via pp7r-overlaps: the 
configuration (a) is lower in energy than (b) due to pp7r- 
interaction. 
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See supplementary material at URL for an illustration of 
the HOMO and LUMO, the spatial variation of the elec- 
trostatic potential in a double chain of passivated arsenics, 
and comparison between the results of geometry optimiza- 
tion for several chain lengths of the (AsH2) n chain. 




FIG. 14. Local charge distribution at a topologically stable pair of an over- and under-coordinated arsenic in a double chain 
of hydrogen-passivated arsenics from Fig. 10 of the main text. The LUMO and HOMO are shown in the top and middle 
panels respectively. The electrostatic potential is shown in the bottom panel, red and blue colors corresponding to positive and 
negative values respectively. 
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FIG. 15. Site dependence of As- As bond lengths in a (AsEh),! chain for several chain lengths is shown, to partially assess the 
sensitivity of the calculation to the chain length, n. (See the main text for calculational details.) For the reader's reference we 
note that the energy gradients in the geometry-optimized chains were about 0.01 kcal/mol/A for even n, and 0.001 kcal/mol/A 
for odd n. 
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